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Perturbative renormalization parameters for heavy quarks* 
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We study the heavy quark mass dependence of the perturbative renormalization factors for the heavy-light 
currents involving Wilson, Clover, and NRQCD heavy quarks. We find that the one-loop Z-factor for the axial- 
vector current depends significantly on the heavy quark mass commonly for all actions we study, while that for 
the vector current has smaller dependence. 
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1. Introduction 

Formulation of the heavy quark on the lat- 
tice, such as NRQCD jjj or the Fermilab ap- 
proach with the Wilson-like relativistic quarks 
requires perturbative matching with the contin- 
uum theory. Since the quark mass in the lattice 
unit cannot be neglected for these theories, the 
perturbative calculation is much more involved 
and have not developed as that for the massless 
case. We report the current status of our effort to 
calculate the one-loop coefficients of the pertur- 
bative Z-factor Zy and Za- It is particularly im- 
portant to calculate the mass dependence of the 
axial-current Z-factor Za, which is required for 
the calculation of the B meson decay constant /b, 
since the one-loop coefficient in the static limit is 
known to be very large 

The combinations of the heavy and light quark 
actions we investigate in this report are listed in 
Table [|. The Wilson and Clover heavy quarks are 
used with the Fermilab group's non-rclativistic 
interpretation with which the relativistic ac- 
tions can be treated in parallel with NRQCD for 
heavy quarks, since the both methods are based 
on the small spatial momentum expansion, or 
equivalently the Aqcd/M expansion. And in 
particular the both methods have the identical 
static limit. 



Table 1 

The combinations of heavy and light quark ac- 



tions. 


Heavy quark 


Light quark 


NRQCD 


Wilson 


Wilson 


Wilson 


Clover 


Clover 



The earlier works for the lattice perturbative 
calculation for NRQCD can be found in Refs. . 
Advanced studies using NRQCD-Clover actions 
are also reported by J. Shigemitsu and A. Ali- 
Khan§. 

2. Operator Matching 

Let us present the operator matching problem 
taking the axial-vector current as an example. 
The matrix element of the time component of 
the axial- vector current Aq for the free heavy and 
light quark external fields having momenta p and 
k have the following form: 



(AsIJn/w/oIp) - t I i (fc)^(C 1 0) +a s Ci iJ )7574 

+ a(C^ 0) + a a d£ ] )l5liil ■ P 

+ aotsC^hj ■ fc7574 

+ higher orders ] h (1) 

where it; is a light quark 4-component spinor, h is 
a heavy quark 2-componcnt spinor, and the coeffi- 
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cients , ■ ■ ■ are functions of the bare quark 
mass o,Mq. The first term gives the leading con- 
tribution surviving in the static limit, and the sec- 
ond and third terms contribute to the Kqcd/^q 
and ah.QCD corrections. 

We have to match the continuum (MS) and the 
lattice theories to give the same matrix elements. 
The matching of the tree level coefficients 
and is almost trivial, and we have calculated 
in this study the one-loop correction , which 
gives the multiplicative Z-factor Za, for all com- 
binations of the heavy and light quark actions. 
The remaining coefficients an d must be 
calculated to achieve higher accuracy. 

3. NRQCD-Wilson 

We use two types of the NRQCD action in- 
cluding terms through 0(Aqcd/M) for the heavy 
quark. The actions are given by 

aHo' 



S = Y,Q(t,x) Q(t,x) 
t,x 

aSH 



1 



2n 



1 - 



2 

2n 



aSH 



Q(t-l,x) 



and 
S -- 



E 



Q(t + l,x) 
aH " 



2n 



type A, (2) 



1 - 



2n 



Q(t,x) 



-Q(t,x)(l-a5H)Q(t,x) , type B, (3) 

where Q is the effective 2-component spinor field. 
H Q = -A^/[2M ] and SH = -gcr ■ B/[2M ] 
are defined as usual. The original 4-component 
spinor field ip is related via Foldy-Wouthuysen- 
Tani transformation through 0(Aqcd/M): 

7 • A N 



i> = i - 



2Mn 



Q \ - 
o 



Rh. 



(4) 



Then the lattice current operators at tree level 
can be written as 

1 



J N = qiTh - 



2M< 



-®r( 7 • A)h 
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Figure 1. One- loop coefficients for the vector and 
axial-vector current renormalization factors for 
the type A action, o: time component of the 
axial vector current; □: time component of the 
vector current; O: spatial component of the vec- 
tor current. Open symbols are obtained with the 
leading term only and filled symbols are with 
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(0) 



(1) 



J, 



(1) 



(5) 



We calculate the multiplicative part of current 
renormalization factors 



Z=l+g 2 (C coat 
where C cont and C latt 



C latt ), 



(G) 



are the one-loop correc- 
tions for the continuum and lattice (axial-)vector 
currents. Figure 1 shows the 1/Mq dependence 



of G 



cont 



C 



latt 



log(aMo)/47T 2 . We have ap- 



plied the tadpole improvement using the average 
plaquette for the mean link variable. The axial- 
vector current in the continuum (MS scheme) is 
defined with the totally anti-commuting 75 . 

We observe that the time-component of the 
axial- vector current and the spatial component of 
the vector current have a large 1/M dependence, 
and in the static limit we reproduce the value 
—0.1346 for the axial- vector current ||. This large 
dependence has an impact on the calculation of 
Jb- For a typical value of the lattice spacing the 
heavy quark mass becomes aMh=1.8-2.5, and the 
difference of Za from the static value could be 
large as ~ 10%. The mass dependence of the 
time component of the vector current is rather 
mild, on the other hand. A similar behavior is 
observed for the type B action. 
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4. Wilson/Clover 

The renormalization factors for the Wilson- 
Wilson and Clover-Clover heavy-light current are 
calculated in a similar manner as the NRQCD- 
Wilson case. The lattice current operator is writ- 
ten at tree level as 



J, 



w/c 



qiTq h + dxqiTij ■ A)q h 



,1 



(0) 



J, 



(i) 



(7) 



where d\ = amo/[2(2 + amo)(l + amo)] with 
amo the bare heavy quark mass appearing in 
the action. In the non-relativistic interpretation 
approach the heavy quark kinetic mass ami is 
given by ami — e ami sinhami/(l+sinhami) and 
ami = log(l + amo). The coefficient d\ is small 
(~ 0.1 or smaller) around the b-quark mass region 
and we neglect in this study. 

The one-loop coefficient for the time compo- 
nent of the vector and axial-vector currents is 
shown in Figure |2| for Wilson- Wilson and Clover- 
Clover as well as NRQCD- Wilson. We apply the 
tadpole improvement by using the critical hop- 
ping parameter to define the mean link variable 
for the relativistic actions. Overall tendencies 
of the mass dependence are similar for all ac- 
tions, except in the lighter heavy quark mass re- 
gion, where the results for the relativistic actions 
tend to have smaller mass dependence while the 
NRQCD- Wilson result still has large slope. The 
mass dependences for the vector current are weak 
compared with that for the axial- vector current. 

5. Summary 

We have investigated the heavy quark mass de- 
pendence of multiplicative part of the renormal- 
ization factors for the heavy-light currents sys- 
tematically using NRQCD, Wilson, and Clover 
actions. We found the large heavy quark mass 
dependence for Za 

The calculation of the heavy-light decay con- 
stant is presented in Ref.[Q for NRQCD and in 
Ref. H for the relativistic actions. 
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Figure 2. The heavy quark mass dependence of 
the renormalization factors for the time compo- 
nent of axial vector currents and vector currents, 
o : NRQCD- Wilson; □ : Wilson- Wilson; O : 
Clover-Clover. Open and filled symbols are for 
the axial vector current and the vector current 
respectively. 
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